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Time-dependent density-functional theory (TDDFT) treats dynamical exchange and correlation (xc) via a 
single-particle potential, 14c(r, t), defined as a nonlocal functional of the density n{r' ,t'). The popular adia- 
batic local-density approximation (ALDA) for T4c(r, t) uses only densities at the same space-time point (r, t). 
To go beyond the ALDA, two local approximations have been proposed based on quantum hydrodynamics and 
elasticity theory: (a) using the current as basic variable (C-TDDFT) [G. Vignale, C. A. Ullrich, and S. Conti, 
Phys. Rev. Lett. 79, 4878 (1997)], (b) working in a co-moving Lagrangian reference frame (L-TDDFT) [I. V. 
Tokatly, Phys. Rev. B 71, 165105 (2005)]. This paper illustrates, compares, and analyzes both non-adiabatic 
theories for simple time-dependent model densities in the linear and nonlinear regime, for a broad range of time 
and frequency scales. C- and L-TDDFT are identical in certain limits, but in general exhibit qualitative and 
quantitative differences in their respective treatment of elastic and dissipative electron dynamics. In situations 
where the electronic density rapidly undergoes large deformations, it is found that non-adiabatic effects can 
become significant, causing the ALDA to break down. 

PACS numbers: 71.10.-w, 71.15.Mb, 71.45.Gm, 73.21.Fg 



I. INTRODUCTION 

Time-dependent density-functional theory (TDDFT) has 
gained considerable popularity as a practical approach to the 
dynamics of many-electron systems.'-^'''* The essential idea 
of TDDFT is to describe N interacting electrons moving in 
an external potential V{r, t) in terms of an auxiliary non- 
interacting system governed by the time-dependent Kohn- 
Sham (TDKS) equation, 
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+ y(r,t) + -(/H(r,t) + Kc(r,i) 



ipair,t) = 



(1) 

(we use atomic units h = e = m = 1 throughout). Here 
and in the following, we consider systems whose ground state 
and dynamical response are everywhere nonmagnetic, and we 
may therefore ignore the spin degree of freedom. Eq. Q thus 
describes the time evolution of doubly occupied Kohn-Sham 
orbitals (pa{j^, t), and the time-dependent density 



N/2 

n(r,i) = 2^|(^„(r,t)|^ 



(2) 



is obtained in principle exactly. In Eq. Q, Vh (r, t) = 
J (Pr'n{r' ,t)/\r — r'| is the time-dependent Hartree poten- 
tial, and Vxc{j^,t) is the exchange-correlation (xc) potential. 
In practice, suitable approximations for X4c(r, t) are required. 
We assume in the following that the system evolves from its 
ground state att = to, although this assumption is not strictly 



Vxc [n] (r, t) has a functional dependence on 



necessary. 
The exact 

n{r,t) that is nonlocal in space and time, i.e., contains in- 
formation about the previous history of the system, includ- 
ing its initial statei^ However, almost all present applications 
of TDDFT employ the adiabatic approximation, ignoring all 
functional dependence of V^c on prior time-dependent densi- 



ties n{r',t'), t' < t. The simplest example is the adiabatic 
local-density approximation (ALDA): 
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where exc [n] is the xc energy density of a homogeneous elec- 
tron gas of density n. The neglect of retardation in ALDA 
implies frequency-independent and real xc kernels in linear 
response.^ This approach has been widely used for calculat- 
ing molecular excitation energies. ^-^ 

One can make the general statement that the adiabatic ap- 
proximation works well for excitations of the many-body sys- 
tem that have a direct counterpart in the Kohn-Sham system, 
such as atomic and molecular single-particle excitationsiS 
On the other hand, for more complicated processes such as 
double or charge-transfer excitations the ALDA often fails 
dramatically.'"'" 

Several attempts to go beyond the ALDA can be found in 
the Iiterature.*'-i2'i_mi5i6j7a8,i9 yign^le and Kohr.'^ showed 

that a non-adiabatic local approximation for exchange and 
correlation requires the time-dependent current \ {y , t) as ba- 
sic variable, rather than the density ri(r, t) (C-TDDFT). This 
formalism was later cast in a physically more transparent 
form using the language of hydrodynamics, '^•^'^ where non- 
adiabatic xc effects appear as viscoelastic stresses in the elec- 
tron liquid. 

To date, C-TDDFT has been applied mainly in frequency- 
dependent linear response. The first application was to 
calculate linewidths of intersubband plasmons in semicon- 
ductor quantum wells. In the absence of disorder and 
phonon scattering, the ALDA gives infinitely sharp plas- 
mon lines. C-TDDFT includes damping due to electronic 
many-body effects, in good agreement with experimental 
linewidths. Van Faassen et alM- calculated static axial 
polarizabilities in molecular chains, which are greatly over- 
estimated with ALDA, and achieved an excellent agreement 
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with ab initio quantum chemical resuhs. Other recent studies 
used C-TDDFT to calculate atomic and molecular excitation 

23.24.25 



energiesii 

Beyond linear response, a wealth of interesting electron dy- 
namics can be explored using TDKS theory. 2'' The C-TDDFT 
formalism has recently been applied to describe linear and 
nonlinear charge-density oscillations in quantum wells in the 
time domain. It was shown that the retardation caused by the 
memory dependence of the xc potential has the striking conse- 
quence of introducing decoherence and energy relaxation.^^-^** 

An alternative non-adiabatic theory has recently been de- 
veloped by one of the authors. '^"^ The idea is to relate the 
local stress in the electron liquid, and thus the xc potential, to 
the dynamics of deformations of fluid elements in the quan- 
tum many-body system. This leads to a formally exact refor- 
mulation of TDDFT from the point of view of an observer in a 
co-moving Lagrangian reference frame (L-TDDFT). Casting 
the theory in terms of Lagrangian coordinates allows one to 
get around the well-known problem of "ultranonlocality" in 
TDDFT, and to derive, in a rigorous fashion, an exact time- 
dependent, non-adiabatic extension of the ground-state LDA. 

In the L-TDDFT formalism the xc potential appears as a lo- 
cal functional of the dynamic deformation tensor At present, 
two limiting forms of this local functional are available. A 
high-frequency, "elastic" form of the non-adiabatic xc po- 
tential was derived in Ref. [iSj. The elastic approximation 
correctly accounts for all complicated nonlinear deformation 
effects, but completely neglects possible xc contributions to 
dissipation. The second available limiting form of the non- 
adiabatic xc potential corresponds to the regime of small de- 
formations. As we will show, in the limit of small deforma- 
tions L-TDDFT formally reduces to C-TDDFT. In this limit 
the xc stress (and thus the xc potential) is proportional to the 
linearized strain tensor which can be considered as a local lin- 
ear functional of the current. In fact, C-TDDFT of Ref. 13 
can be viewed as a viscoelastic linear Hooke's law in the non- 
linear quantum continuum mechanics defined by the general 
formulation of L-TDDFT. In contrast to the purely elastic ap- 
proximation constructed in Ref. C-TDDFT contains all 
dissipation effects, but is formally restricted to infinitesimally 
small deformations. Of course, the formal asymptotic criteria 
and the practical regimes of applicability of any approxima- 
tion can be quite different. Therefore a more detailed analysis 
based on explicit numerical examples is required to assess the 
validity of the two currently available non-adiabatic xc func- 
tionals, and to analyze their relation to ALDA. 

Thus, the purpose of this paper is to illustrate, compare, 
and analyze C-TDDFT and the elastic approximation to L- 
TDDFT for simple, quasi-one-dimensional model systems, in 
order to show the differences and common grounds of both 
approaches. We explore the performance of the two xc poten- 
tials for two kinds of analytically given time-dependent model 
densities, representing charge-density oscillations in the form 
of collective sloshing and breathing modes of varying ampli- 
tudes and frequencies. This will allow us to simulate electron 
dynamics over a wide range of time and frequency scales, 
from the linear to the nonlinear regime. In our analysis, we 
will focus on a detailed comparison of the time-dependent xc 



potentials, as well as on the instantaneous and time-averaged 
power absorption associated with the charge-density oscilla- 
tions. This will give us insight into the inner workings of C- 
and L-TDDFT in different dynamical regimes, and their rela- 
tion to the ALDA. In particular, we will discuss and clarify the 
meaning of "non-adiabatic", the cross-over from the linear to 
the nonlinear domain, and the competition and coexistence of 
elastic and dissipative xc effects. 

It turns out that, for small-amplitude deformations, C- and 
L-TDDFT agree in the limit of short time scales or, equiv- 
alently, in the high-frequency regime. In general, the two 
theories show some differences in their treatment of elastic 
and dissipative effects in the inhomogeneous electron liquid. 
However, the size of these differences strongly depends on the 
type of collective mode, and the associated charge-density de- 
formations. If the deformations are large and occur on short 
time (high frequency) scales, both non-adiabatic theories give 
a clear indication of a failure of the ALDA. 

This paper is organized as follows. In section HIl we sum- 
marize the essential formal framework of C- and L-TDDFT. 
Section|lIl]shows how to construct simple analytic model den- 
sities in the Lagrangian and associated laboratory reference 
frame. Section HVl presents detailed numerical results and dis- 
cussion, with a separate treatment of the linear and nonlinear 
regime. We give our conclusions in SectionfVl 



II. NON-ADIABATIC THEORIES IN TDDFT 



A. C-TDDFT 



1. Linear-response regime 



Starting point of C-TDDFT is the linear current-density re- 
sponse ji(r,cj) to an external, frequency-dependent vector 
potential Acxti(i", t^). 



ji,,{r,uj) = d^r'xKS,ij{r,r',uj)[AcM,j{r',uj) 



(4) 



where i,j denote Cartesian coordinates, and XKS,ij is the non- 
interacting, Kohn-Sham current-current response tensor. Here 
and in the following, we use the Einstein convention for the 
summation over repeated indices. The Hartree vector poten- 
tial is given by 



AHij(r,Lj) 



, V'-ji(r',c^) 
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The simplest approximation for the linearized xc vector po- 
tential Axel is the ALDA, which is defined as 
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■y'dV/-°*(r,r')V'-ji(r',u;), (6) 



where 
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rt=no(r) 
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is the frequency-independent ALDA xc kernel [tiq (r) is the 
ground-state density]. In contrast with the xc scalar poten- 
tial, the xc vector potential admits a frequency-dependent lo- 
cal approximation. The resulting expression can be written as 
follows; '^i^^i'^ 



iLunQ{r) 



Vfccrxcjfc(r,tj) . (8) 



Here, c denotes the speed of light, and cTxcjk is the xc vis- 
coelastic stress tensor: 



xc,j/c 



VjWi,fc + VfeWi J - -V • vi Sjk 
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CxcV • Vi Sjk 



(9) 



where Vi(r,cj) = ji(r,CL')/no(r) is the velocity field associ- 
ated with the current response, and rj^c and (^c are complex 
viscosity coefficients defined as 



(10) 



Cxc("-,'^) 
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f^^{n,u!) and f'^^{n,uo) are frequency-dependent xc kernels 
for the homogeneous electron gas, which can be found in var- 
ious parametrizations in the literature . ^i^^i-^" In Eq. -q^c 
and ^xc are both evaluated at the local no(r). 



2. Nonlinear regime 

The generalizationi^iS of C-TDDFT into the nonlinear 
regime and the time domain requires solving the following 
TDKS equation: 



dt 



■ A(r,<) + - Axc(r,t) 



+ V{v, t)+Vn{v,t) 



(12) 



Notice that the Hartree term can be expressed as a scalar po- 
tential, and we are free to admit external scalar as well as vec- 
tor potentials. 

As explained in Ref. the form of the nonlinear xc vec- 
tor potential is dictated by a number of general requirements, 
such as Newton's third law (xc force density follows from a 
symmetric stress tensor), and the proper limit in the linear 
regime, which was discussed in the subsection above. A for- 
mally exact, general expression for Axe resulting from these 
requirements will be presented within the Lagrangian frame- 
work in Section HTbI However, a straightforward expression 
for a nonlinear, nonadiabatic xc vector potential, valid up to 
second order in the spatial derivatives, follows almost imme- 
diately: 



1 dA^,^,{v,t) 
c dt 



n{r,t) 



, (13) 



where the viscoelastic stress tensor CTxc now contains the time- 
dependent velocity field v(r, t) — j(r, t)/n{r, t): 



^xc,ij (j^ : ^) 



dt' \7ir, t,t') 



\/^Vj{r,t')+VjVi{r,t') 



+ C(r,i,t')V-v(r,t')% M14) 



-V.v(r,t')% 



The viscosity coefficients in Eq. ( I14t are the Fourier trans- 
forms of ilO\ and il 1> : 



ri[r,t,t') = / ^ v{n,uj)e' 



j{t~t') 



(15) 
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and similar for The apparent ambiguity in Eq. ( I15> whether 
the density should be evaluated at t or t' is resolved by noting 
that the difference involves higher gradient corrections. We 
emphasize again that the simple form of Eq. jl3t is justified 
by the basic assumption that the gradients of the density and 
velocity are small; the velocity itself, on the other hand, need 
not be small. These points will be elaborated in more detail 
in Section HTBl where we will explain how the approximate 
expression ( I13> is obtained from the Lagrangian framework 
in the appropriate limit. 

In the following, we consider model systems where all spa- 
tial dependence is along the x direction only. One can then 
transform the xc vector potential, Eq. ( I13L into a scalar one: 
V^cix,t) = V^^^^{x,t) + V^{x,t) (ALDAh-M), with the 
memory part given by 



V^{x,t) 



dx' 



n{x' , t) 



{x',t) 
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Assuming that the system has been in the ground state (with 
zero velocity field) for t < 0, the xx component of the xc 
stress tensor becomes 



(17) 



<Jxc,xx{x',t)= / Y{n{x',t),t-t')V,,v,,{x',t')dt' 
Jo 

where the memory kernel Y is given by 



Y{n, t-t') = - ri{n, t ~ t') + ({n, t - t') . (18) 
o 

With the help of the Kramers-Kronig relations for f^^ we can 
express the memory kernel as follows: 



Y{n,t-t') = -^x 



n 
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5/i;(u;)cosLo(i-t'), (19) 



with the xc shear modulus of the electron liquid. 
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The short-time behavior of Y{n, t — t') is of particular inter- 
est, since it governs the high-frequency dynamics. The limit 
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FIG. 1: Scaled memory Kernel Y{n, t - t') for = 1, 2, 3, 4, 5, 
using the GK (full lines) and QV (dashed lines) parametrizations. 
Inset: short-time limit Y{n, 0) [Eq. <22H versus fs- 



the convective motion of the electron fluid (the particles at 
a given point of space retain the memory of their previous 
positions)."* The key idea of L-TDDFT is to eliminate the 
above source of nonlocality by reformulating the theory in 
the Lagrangian frame, i.e., in a local reference frame moving 
with the fluid. Since the convective motion in the Lagrangian 
frame is absent, a spatially local description of xc effects be- 
comes possible. This possibility represents the most important 
outcome of L-TDDFT: it allows one to derive an exact nona- 
diabatic extension of LDA into the dynamic regime. 

The general formulation of L-TDDFT starts with the ex- 
act relation of the xc force to the xc stress tensor P^cij'^^ 
By definition the xc potential Axe ensures that the physical 
density and current are reproduced by an auxiliary system of 
noninteracting KS particles. This means that Axe should pro- 
duce an effective xc Lorentz force that exactly compensates 
for a difference of local stress forces in the real interacting 
system and in the noninteracting KS system. Accordingly, the 
xc vector potential should satisfy the following equation: 



Y{n,0) can be expressed analytically using the Kramers- 
Kronig relation 
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where the high-frequency limit is known via the third- 
moment sum rulesi^ The result is 
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Y{n,0) = -— exc + — -r- 
6 o an 



It is also straightforward to see that ^{n,0) = 0, i.e., the 
memory kernel starts with zero slope. 

Figure n shows the memory kernel Y{n,t — t') evalu- 
ated with the Gross-Kohn (GK)^ and Qian-Vignale (QV)^" 
parametrizations for f^^ [lo), and scaled by the short-time limit 
y(n,0) (see inset). Here, Tp — 2ti/lOp is the characteristic 
time scale associated with the plasma frequency Up = ^/Airn. 
As noted earlier, yGK yQV ^.g j^jentical att = t' and 
have a similar short-time behavior. Y'^^ goes exponentially 
to zero for large t — t', and the scaled GK memory kernels are 
nearly identical for all times over a wide range of r^. On the 
other hand, all Y'~'^ pass through a minimum around 0.22Tp 
and then approach the finite limit 4/Xxc/3, where /Xxc 
for large r^i^ It is interesting to notice that both QV and GK 
memory kernels reach their long-time asymptotic limits very 
rapidly, within about half a plasma cycle Tp. 



B. L-TDDFT 

1. General formulation 

The main problem in constructing nonadiabatic approxima- 
tions for xc potentials is an inherent nonlocality of the time- 
dependent theory. Physically, this nonlocality is related to 



at 



where Px, 



xc.ij — ^ ~ T]^^ is the difference of the full stress 
tensor for the interacting system, Pij, and the kinetic stress 
tensor for KS system, Tf^^. Equation ( I23> serves as a basic 
definition of Axe, which automatically accounts for the zero 
force and zero torque conditions.'^ 

We note that the xc stress tensor Pxc,ij is to be distinguished 
from the earlier introduced xc stress tensor a^cij- The con- 
nection between the two will be explained in detail in subsec- 
(22) tion lll B 2l The main difference lies in the fact that P^c.ij, for- 
mally exactly and to all orders in the inhomogeneity, accounts 
for all dynamical xc effects, whereas the ALDA has been sep- 



arated out in the definition of cr 



Furthermore, a^c,ij is 



valid only for small deformations of the electron liquid (in a 
sense to be defined below). 

All the rest of L-TDDFT can be viewed as a calculation 
of the dynamic xc stress tensor P^c,ij, which enters the defi- 
nition of Axe, Eq. i23\ . by reformulating the problem in the 
co-moving Lagrangian frame. The transformation to the La- 
grangian frame corresponds to a nonlinear transformation of 
coordinates, r = r(^, t), where r(^, t) is the trajectory of an 
infinitesimal fluid element that evolves from the point ^. For- 
mally the function r{^,t) is defined by the following initial 
value problem: 



at 



(24) 



The initial positions, ^, of the fluid elements play the role of 
spatial coordinates in the co-moving frame. The above trans- 
formation from the old coordinates r to the new coordinates 
^ induces a change of metric. The metric tensor in the La- 
grangian ^-space coincides with the Green's deformation ten- 
sor. 



arkitt) aMtt) 
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(25) 



5 



which is a common characteristic of deformations in the La- 
grangian formulation of continuum mechanics. A complete 
reformulation of TDDFT in the co-moving frame shows that 
a local description of xc effects is possible if one uses t) 
as a basic variable: the stress tensor in the Lagrangian frame 
can be consistently considered as a spatially local functional 
of the Green's deformation tensor. 

For most practical applications we need the xc potential in 
the laboratory frame, Eq. i23\ . Transforming the stress tensor 
functional from the co-moving to the laboratory frame, we 
find the required tensor P^c,ij entering Eq. ( I23> . The locality 
in the Lagrangian frame translates to a local dependence of 
Pxc,ij on the Cauchy's deformation tensor. 



dri 



(26) 



where the function ^(r, t) is obtained by inverting the trajec- 
tory equation, r = r{^,t). The Cauchy's tensor gij{r,t) is 
commonly used to describe deformations in the Eulerian for- 
mulation of continuum mechanics. From a physical point of 
view, the functional dependence of the stress tensor on gij thus 
emerges quite naturally. 

Substitution of Pxc,ij[gij] into Eq. (I23> yields the xc vec- 
tor potential A^c as a functional of the Cauchy's deformation 
tensor. If gradients of the deformation tensor are small, the 
functional Pxc,ij[9ij] is local in space, but in general it can 
be nonlocal in time. The time-nonlocality of the stress tensor 
functional is related to the memory loss due to dissipation. In 
the absence of dissipation (e. g. in the exchange-only approx- 
imation) the stress tensor becomes a simple function of gij, 
which corresponds to an infinitely long memory. 

Equations i24i and (I26> determine the Cauchy's deforma- 
tion tensor as a functional of velocity v(r, i). An alterna- 
tive, and practically more convenient way to compute for a 
given velocity is to solve the equation of motion that governs 
the dynamics of (r, t) directly in the laboratory frame: 



dt 



dg. 



dv, 



dvk 



^fc — = - g^rakj -—9^k, g^J (r , O) = (5y 



drk 



(27) 

An important property of the deformation tensor is that it al- 
lows us to relate the time-dependent density n(r, t) to the ini- 
tial density distribution, no(r): 



i{r,t) = ^/g{r,t)noi^{r,t)), 



(28) 



where g{r, t) is the determinant of gij (r, t). 

Equation (I27> or, equivalently, Eqs. (I24> and i26\ . show that 
in general the deformation tensor is a strongly nonlocal (both 
in space and in time) functional of the velocity. Therefore, in 
spite of the fact that the xc stress tensor and consequently the 
xc vector potential are local functionals of cjij, they are nonlo- 
cal in terms of velocity or any other variable. This makes the 
choice of (r, t) as a basic variable much more preferable. 

To obtain an explicit construction of the local functional 
Pxc,ij[gij], a solution of a homogeneous time-dependent 
many-body problem in the Lagrangian frame is required (see 
Sec. V in Ref..l8.). A complete solution of this problem seems 



to be impossible, at least at the current level of knowledge. 
However, there are two practically important, exactly solvable 
special cases, which are described below. 



2. Small deformation approximation: Recovery of C-TDDFT 

The abovementioned many-body problem in the homoge- 
neously deformed Lagrangian ^-space can be solved pertur- 
batively if the deformation tensor gij only slightly deviates 



from the unit tensor 8, 



V 



gy(r,t) = S,j + 6g,j{r,t) 



(29) 



Introducing the displacement vector, u(r, t) = r — ^(r, t), and 
using Eq. ( I26t . we find that small Sgij corresponds to small 
gradients of the displacement: 



Sgij{r,t) = - 



dui 
dri 



dri 



(30) 



Clearly, small gradients of u(r, t) imply that the velocity gra- 
dients are also small, since to lowest order in WiUj Eq. ( I24t re- 
duces to the relation 9tu(r, t) = v(r, t). Obviously the small- 
ness of deformations does not mean that the displacement or 
the velocity themselves are small (i.e., the system can be far 
beyond the linear response regime). A well known example is 
the rigid motion of a many-body system in a harmonic poten- 
tial, where gij = Sij but the displacement can be arbitrarily 
large. 

The stress tensor functional for small displacement vectors 
If— » ^ 
was derived in Ref. 18. Extension of this derivation to the gen- 
eral regime of small deformations, i.e., to the regime of small 
displacement gradients, is straightforward. The resulting xc 
stress tensor takes the following form: 



(31) 



where P^^^{n) is the xc pressure of a homogeneous elec- 
tron gas, and SPij is a nonadiabatic correction, which is linear 
in Sgif. 

t 

SP,,{r,t) = J dt'I^^K^,{n{r,t),t-t')Sgkk{r,t') 



(32) 



The kernels fi^dn, t — t') and K^dn, t — t') in Eq. \3>2\ have 
the meaning of nonadiabatic shear and bulk moduli, respec- 
tively [the adiabatic part of the bulk modulus is included in the 
ALDA pressure term in Eq. \3\V\. The corresponding Fourier 
transforms of the elastic moduli, fixcin, lli) and Kxc{n, w), are 
related to the complex viscosity coefficients of Eqs. JIOI I and 
(II 1> . rixc{n,u!) and (xcinjio), as follows: 



lJ-xc{^) = -iw?7xc(t^), Kxc{uj) = -iwCxc(w) 



(33) 



Using Eq. ( I33> and the relation dtu ~ v, we find that the 
nonadiabatic stress tensor SPij, Eq. (l32t . is identical, up to a 
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sign,"" to the tensor a^c.ij of Eq. ( I14> [i. e. 5Pij = —cr^cAj}- 
In addition, in the Hmit of small displacement/velocity gra- 
dients the spatial derivatives of Axe on the left hand side of 
Eq. \22>\ are negligible. Thus, in the regime of small de- 
formations we recover the complete nonlinear form of C- 
TDDFTji^ Eqs. Q and (O. 

The imaginary parts of the complex elastic moduli, K^c (^^) 
and /ixc(ijj), are responsible for the dissipation (viscous) ef- 
fects. For the high-frequency/short-time dynamics these ef- 
fects become irrelevant. As a result, the high-frequency limit 
of the nonadiabatic stress tensor of Eq. ( I32t becomes com- 
pletely local and purely elastic: 



-^Sgkkir, t) 



(34) 



where (n) and /i^ (n) are the high-frequency limits of the 
bulk and shear moduli, respectively. 

The structure of the small deformation approximation, 
Eqs. (13 U and i32\ . clearly demonstrates that in this regime 
the nonadiabatic contribution appears as a small, linear in Sgij 
correction to the adiabatic dynamics. If the process is strongly 
nonadiabatic, the deformations can not be considered small. 
In fact, the deviation of the deformation tensor from Sij can 
serve as a general measure of nonadiabaticity. 



3. Nonlinear elastic approximation to L-TDDFT 

It is very difficult to account both for the full nonlinear de- 
pendence on gij, and for the dissipation. C-TDDFT includes 



all xc dissipation effects on a level linear in 



On the 



contrary, if we neglect the dissipation effects, a closed non- 
linear local approximation for the stress tensor can be rigor- 
ously derived.'^ The reason is that the homogeneous many- 
body problem, which has been formulated in Ref. 18, admits a 
simple complete solution in the regime of fast dynamics when 
the dissipation is irrelevant. In this case the xc stress tensor 
Pxc,ij (r, t) becomes a function of the time-dependent density 
n(r, t) and the Cauchy's deformation tensor gij (r, t): 

Pxc,ij = -^gijVdE^fn (^'^^ + Lijigki) Spot (^-^ 

(35) 

where E^f^{n) and £'pot("-) are the xc kinetic energy and the 
potential energy per unit volume of the homogeneous electron 
gas. The function Li^igki) in Eq. (I35> is explicitly defined in 
Appendix A of Ref. 1 8 . In the limit of small deformations the 
nonlinear elastic approximation of Eq. (I35> can be expanded 
around gjj — 6ij and reduces to the linearized form defined 
by Eqs. OH and (I34> . In other words one recovers the high- 
frequency limit of C-TDDFT. 

We conclude this section with the explicit formulation of 
the nonlinear elastic approximation for a one-dimensional 
motion. If all spatial variations are along the x-axis only, the 
deformation tensor takes a diagonal form with g^z = gyy = 1, 



and gxx = g{x, t). The xc effects can then be described by an 
xc scalar potential that is related to the xc stress tensor as fol- 
lows: 



dx' d 



n{x' ^ t) dx' 



',t)rg{x\t)). (36) 



Equation ( I35> for the xx-component of the xc stress tensor 
reduces to the form 



Pxc,..(n,5) = \f'^El^^ + L{g)E, 



pot 



n \ 

71) 



where the factor L{g) is given by 



m = 



5-1 



1 - 



arctan ^/g^T 



(37) 



(38) 



Finally, Eq. ( I27> . which relates the deformation g{x, t) to the 
velocity v{x, t), simplifies as follows: 



dg_ ^ _^^dg_ _ ^dv__ 
dt dx dx^^ 



7(x,0) = 1. 



(39) 



It is worth mentioning that L{g 1) = 1/3 in the limit of 
zero deformation. Eq. (I37> then approaches the standard virial 
expression for the xc pressure, and V^^{x^t) reduces to the 
ALDA xc potential. We define the post-ALDA contribution 
of the nonlinear elastic L-TDDFT xc potential as 



V^,{x,t)^V^,{x,t)-V^^''^{x,t) 



(40) 



To shorten the notation for the rest of this paper: from now 
on, whenever we refer to L-TDDFT, we actually mean the 
nonlinear elastic approximation to L-TDDFT, defined in this 
section. 



III. CONSTRUCTION OF ONE-DIMENSIONAL 
ANALYTICAL EXAMPLES 

In one dimension, the trajectory of a fluid element with La- 
grangian coordinate is determined by the following equa- 
tion: 



dt 



(41) 



X is the position, at time t, of a fluid element that started out 
at ^, and v is its velocity. In general, this is a complicated 
nonlinear differential equation, with formal solution 



(42) 



In other words: if we know the velocity, at time t, of a fluid 
element that started out at ^, we can determine its trajectory by 
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direct integration. From this, we can then determine the time- 
dependent density in the laboratory frame: we first invert (142 > 
to obtain £^{x, t), then compute the deformation as 



and finally obtain 



(43) 



(44) 



In practice, of course, this procedure is not very helpful, since 
the functional form of v{^,t) is unknown. However, we can 
use it to construct simple analytic examples, as follows. 



A. Sloshing mode 

We assume that the system is confined within hard walls, 
-L/2 < (x,^) < L/2, with initial density 



naiO = — cos^ ( Y 



(45) 



where N is the number of electrons per unit area (sheet den- 
sity) in the i/ — z plane. We assume a simple quadratic form 
of the velocity field: 




FIG. 2: Snapshots of density n{x,t), in units of N/L, velocity 
v{x, t), in units of Lcj, and deformation t) for the sloshing mode 
in the laboratory frame, taken at times t = 0, r/4, T/2, 3r/4. 



t) = Alo ( ^ ^ ) coswt . 



Equation J42> is then easily integrated: 



(46) 



t) =^ + A [ ^ - J- \ sintot . 



(47) 



The next step is to invert Eq. i47l to determine the trajecto- 
ries of the fluid elements, which requires solving a quadratic 
equation, with the result 



L 



2 A sin ujt 



_ 9 4Aa; 
1 + A-' sm uit — smcut 



L 



(48) 

which properly reduces to ^ = a; for A 0. The range of 
allowed amplitudes is |^| < 1, which is dictated by the con- 
straint that no fluid element can cross the hard- wall boundaries 
at ±L/2. 

We can now calculate the deformation using Eq. (I43> : 



{x,t) = [1 + A^ shr tut 



4Ax 



sin Lot 



(49) 



The time-dependent density of the sloshing mode in the labo- 
ratory frame, n{x, t), then follows from Eq. i44\ . using (I48> 
and(|49}. 



3 




-0.75 -0.5 -0.25 0.25 0.5 0.75 

x/L 

FIG. 3: Same as Fig.|2| for the breathing mode. 
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B. Breathing mode 

To simulate a breathing mode, we assume a linear velocity 
distribution of the fluid elements: 

v{^,t) ^ Alo^ cos Lot . (50) 

According to Eq. J42t . this gives the following trajectory: 

x{^,t)=^{l + Asmujt), \A\<1. (51) 
This is easily inverted: 

^(x,^)= ^^f. (52) 
1 + A sm uit 

and the resulting deformation is 



(1 + A sin Loty 



(53) 



We choose the same initial density distribution no(^), Eq. 
(I45> . as for the sloshing mode, and the resulting time- 
dependent density of the breathing mode is 



n(x,t) = ——cos — 



{l + Asinujt)-\ (54) 



where |a::| < {L/2){1 + Asmojt). 

Figs.|2land|3lshow snapshots of n{x, t), v{x, t), and g{x, t) 
for the sloshing and the breathing mode, taken at time t — 
0, T/4:, T/2, 3T/4, where T = 2tt/uj. The amplitude is A ^ 
0.5 in both cases, length is measured in units of L, and density 
and velocity are plotted in units of N/ L and Luj, respectively. 
The deformation g{x,t)is maximal at the turning points of the 
osciflations (at t = T/4 and 3T/4). We find that the breathing 
mode features large deformations everywhere, i.e., g deviates 
strongly from one. The sloshing mode, on the other hand, is 
strongly deformed only towards the edges where the density is 
small. We will see below how this affects the non-adiabaticity 
of the xc potential of the two modes. 



IV. RESULTS AND DISCUSSION 

In the following numerical examples, we choose a system 
where, in atomic units, we have sheet density = 1 a.u. and 
quantum well width L — 10 a.u. We shall measure frequen- 
cies in units of the average plasmon frequency of the system, 
given by 



UJr, 



dx ujp{n{x)) 



(55) 



For our initial density distribution nQ{x), Eq. (I45> . we find 

Cjp ~ ^32A/7ri, which for the above values of N and L 
comes out as cop — 1.009 a.u. 




FIG. 4: (Color online) Snapshots of ALDA xc potential fluctuations 
and xc memory contributions during one cycle of the sloshing mode 
(Fig. |2j, with amplitude A = 0.005. Blue solid lines: V^^°^, 
scaled by 0.1 (independent of frequency). Red solid lines: V^'^, in 
GK parametrization, for low frequency, cijl ~ 0.001 ujp, crossover 
frequency, oJc = i-7 CUp, and high frequency, luh = 1000 tjp. Green 
dashed lines: at 0.1, 0.5, 2, 10 times ujc- 



A. Linear regime 

1 . Time-dependent xc potentials 

We first consider the linear regime of small density fluc- 
tuations. For the time-dependent densities associated with the 
breathing mode and the sloshing mode, we calculate and com- 
pare three different time-dependent xc potentials: V^^{x,t) 
from C-TDDFT [Eq. (O], V^^{x,t) from L-TDDFT [Eq. 

and the ALDA potential fluctuations V^^^^{x,t) = 
V^^^^^{n{x,t))-V^^J'^{noix)). 

Figures EJQ each show four snapshots of V^{x,t) and 
V^^^^ during one cycle of the sloshing/breathing modes, at 
t = 0,T and T/2 (density passes through equilibrium, with 
maximal velocity) and t — T/A and 3T/4 (density at turn- 
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FIG. 5: (Color online) Same as Fig. El but V^^ calculated 
parametrization. The crossover frequency is loc = 2.15 uop. 
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FIG. 6: (Color online) Same as Fig. |4| but for the breathing mode 
(Fig. EJ. V^^^ is not scaled. The crossover frequency is loc = 



2.22 u)^ 



ing points), as illustrated in Figs. I2I3I V^{x,t) has been 
calculated using the GK and QV parametrizations as input 
(Figs. 1416 1 and Figs. I5I7I respectively). All density oscilla- 
tions have the same small amplitude A = 0.005, but differ- 
ent frequencies. We highlight the low- and high-frequency 
limit, wl = 0.001 (Dp and wh = 1000 (Dp, respectively, and 
an intermediate crossover frequency uc, to be defined below, 
which varies between 1.7 (Dp and 2.65 Dp for the four cases 
considered. To illustrate the strong dependence on frequency 
in the crossover regime, we also show at 0.1,0.5,2, and 
10 times luq- In the high-frequency limit wh, y^{x,t) is 
identical to V^{x, t) in all cases, i.e., the dynamics is purely 
elastic. 

The results in Figs. |4j0reveal the following features: 

(i) The high-frequency limit ofV^{x, t) [i.e., V^{x, t)] is 
phase-shifted by vr with respect to V^^^^{x,t) in all cases 
considered. This is to be expected for a purely elastic poten- 
tial: it reaches its maximum at the instant of largest displace- 
ment from equilibrium. 

(ii) At the low-frequency limit wl, the behavior of de- 



pends on the parametrization used to calculate the memory 
kernel: in GK, V^'^ vanishes, whereas in QV, again be- 
comes purely elastic. This reflects the different long-time be- 
haviors of the memory kernel: in GK, it decreases exponen- 
tially, whereas in QV it approaches a finite constant (see Fig. 

(iii) At intermediate frequencies, the phase shift between 
yALBA ^jjj yaries between tt and n/2, where tt/2 indi- 
cates a purely dissipative potential (see below). 

(iv) The average strength of the forces associated with 
grows with frequency, and becomes comparable to the 

ALDA fluctuating forces in the high-frequency limit (notice 
that V^^°^ is scaled by 0.1 in Figs. |4|and|5}. This clearly 
shows that non-adiabatic effects can become non-negligible in 
practice. We will say more about this below when we discuss 
the nonlinear regime. 



10 



0.001 
-r 0.0005 

1 
(2 -0.0005 
-0.001 

0.0005 

1 
c 

B 

S. -0.0005 
-0.001 

0.0005 

1 
c2 -0.0005 

-0.001 
^ 0.0005 
3 
<£ -0.0005 

-0.001 




ALDA 



L3 




ALDA 



L,H 




2.5 





X (a.u.) 

FIG. 7: (Color online) Same as Fig. |6| but V^^ calculated in QV 
parametrization. The crossover frequency is cue = 2.65 Cup 



2. Power and dissipation 

For a more quantitative analysis, it is useful to consider the 
power associated with the charge-density oscillations of the 
two types of modes. We define the power in the usual way as 
force density times velocity: 



r{t) 



dxv{x,t)n{x,t) ^V^{x,t) 



(56) 



Figure |S] shows V^t), scaled by utA^, during one cycle of 
the sloshing and breathing modes (calculated using the GK 
parametrization). In the low-frequency limit, where the cur- 
rents are vanishing, the power tends to zero, but as the fre- 
quency increases, a periodic input/output of power is observed 
to take place during a cycle. It can be clearly seen that, on 
average, V is more negative than positive for intermediate fre- 
quencies, which indicates net power dissipation. In the high- 
frequency, elastic limit, V{t) has sizable amplitudes, but aver- 
ages to zero during a cycle. 

These findings are summarized in Figures|9]and[l0| We plot 
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FIG. 8: Power V{t) [Eq. ^] calculated using GK, associated with 
the sloshing and the breathing mode, during one cycle of the charge- 
density oscillation (T = 2-n /uj), for different frequencies. L,C,H 
indicate the low-frequency (purely dissipative in GK), crossover, and 
high-frequency regime (purely elastic). 



~V, the absolute value of the time-average of 'P(t), scaled by 
LoA^, which represents the net power absorption per cycle, for 
the sloshing and the breathing modes. Notice that —V in all 
cases has a pronounced enhancement for frequencies of the 
order of ujp. We define the crossover frequencies ujc as those 
frequencies where the maxima of power absorption occur. 

While the two modes behave qualitatively very similarly as 
far as their power absorption is concerned, one observes in 
Figs. ISHTOI that the absorbed power of the breathing mode 
is about an order of magnitude higher than for the slosh- 
ing mode, for the same value of the amplitude, A = 0.005. 
This hardly comes as a surprise: our sloshing mode can be 
viewed as a cousin of Kohn's mode,^^-^^ and the Harmonic 
Potential Theorem.^^ tells us that charge-density oscillations in 
parabolic quantum wells are undamped. The breathing mode, 
on the other hand, bears no resemblance at all to Kohn's mode. 
The sloshing mode has a smooth hydrodynamic flow with rel- 
atively little internal compression except at the turning points. 
By contrast, the defining feature of the breathing mode is the 
periodic compression and rarefaction of its density profile, 
i.e., a very "un-hydrodynamic" behavior. 

It turns out that there is a direct relation between the av- 
erage power absorption and the imaginary part of the xc 
kernel, ^/^^{uj). This is clear from Eq. (I19> . which ex- 
presses the memory kernel Y{n,t — t') via Fourier trans- 
form of Q /^^{lli) / u! . The full lines in Figs. l9l and [Tol plot 
^T^/xcC"^)' where 7 is a constant scaling factor. For a best 
fit, 7si = 0.00525 and rf = 1.41 for the sloshing mode, and 
7br — 0.071 and r^'' = 1.18 for the breathing mode (same 
for GK and QV). Notice that the equilibrium density no{x) 
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FIG. 9: Full symbols: net rate of power absorption, —V/loA'^, dur- 
ing one cycle of the sloshing and breathing modes, versus frequency 
(in GK). Full hne: -79/4(t^), for rs = 1.41 (sloshing) and 1.18 
(breathing). Open symbols: Phase lag (in units of n) between 
and ALDA. The crossover frequency is defined by a maximum in 
power absorption and a phase lag 3n/4. The dashed lines are a guide 
to the eye. 



has a value of — 1.06 in the center. It thus emerges that 
the dissipation is dominated by different regions of the den- 
sity distribution for the two modes: around x = ±0.27i for 
the sloshing mode, i.e., halfway between center and walls, and 
X = ±0.18L for the breathing mode, i.e., much closer to the 
center. 

The phase lag of the modes as a function of frequency was 
akeady discussed in the context of Figs. I4I7I However, in 
Figs. |9]and^|one can see the difference between GK and QV 
most dramatically. We consider the phase lag between the in- 
stantaneous power absorption, and the ALDA potential 
fluctuations (which are in phase with the charge-density oscil- 
lations). Fig. |9] shows that in GK one has a transition from 
a purely dissipative behavior in the low-frequency limit, with 
phase difference 7r/2, via a crossover region of mixed dissi- 
pative/elastic behavior, to the high-frequency, purely elastic 
regime with phase lag tt. The maximal power absorption oc- 
curs for a phase lag of 37r/4. On the other hand, in the QV 
parametrization the low-frequency regime is also purely elas- 
tic, i.e. both the low- and high-frequency limits have phase lag 
TT, and dissipative contributions come in only at intermediate 
frequencies. The resulting power loss as a function of fre- 
quency is thus a bit higher and more narrowly peaked about 
u!c in QV than in GK. 

As mentioned earlier, the qualitative differences between 
GK and QV have their origin in the different long-time be- 
havior of the memory kernels, see Fig. [^and the discussion in 
Section lTl A2I The overall result is that there is a broad range 
of frequencies, between about 0.1 to 10 times the characteris- 
tic average plasma frequency u)p, where the system exhibits a 
mixed elastic/dissipative behavior, which can lead to substan- 
tial dissipation. 
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FIG. 10: Same as Fig. |9| but in QV parametrization. Both the low- 
and high-frequency limits are purely elastic. 



B. Nonlinear regime 

1. C- versus L-TDDFT in the high-frequency limit 

In the following, we will extend our numerical studies of 
C- and L-TDDFT into the nonlinear regime. We begin by 
directing our attention to the high-frequency, purely elastic 
region, since our approximate version of L-TDDFT becomes 
exact in that limit and can thus be used as a benchmark to 
assess the accuracy of nonlinear C-TDDFT. 

Fig. II II shows the scaled power 'P{t)/LuA'^ [Eq. ( I56H for 
u! = W^LUp during one cycle of the sloshing and breath- 
ing modes, comparing small amplitudes {A = 0.005) and 
large amplitudes (A = 0.5 and 0.75). In the small-amplitude 
regime, 'P{t)/u!A'^ has a sinusoidal shape for both modes, in- 
dependent of the amplitude A of the oscillations (see also Fig. 
|8}, but for large amplitudes, nonlinear deviations occur. In 
the sloshing mode, tends to a more sawtooth-like shape. 
In the breathing mode, V{t) becomes suppressed when the 
charge density spreads out in the first half of the cycle, and 
dramatically enhanced when the density gets squeezed in the 
second half of the cycle. As we will see below, these strong 
deformations give rise to sizable non-adiabatic effects. 

For small amplitude, A = 0.005, C- and L-TDDFT are 
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FIG. 11: Power [Eq. |56}] calculated using GK, over one cycle 
of the sloshing and the breathing modes, for uj = W^ujp. Full and 
dashed lines: C- and L-TDDFT for amplitudes A = 0.005, 0.5, and 
0.75. At A = 0.005, C- and L-TDDFT are identical (see Fig.|8}. 




FIG. 12: Percent deviation of C-TDDFT from the (exact) L-TDDFT 
result for the time-averaged absolute value of the scaled power, {t2 — 
t2)^^ J^*^ dt\V{t)\/u}A'^. Top panel: sloshing mode, averaged over 
one cycle. Bottom panel: breathing mode, averaged over the first 
(full symbols) and the second half-cycle (open symbols). 



identical, but for larger amplitudes, A = 0.5 and 0.75, some 
differences develop. However, the magnitude of these devi- 
ations strongly depends on the type of mode. For the slosh- 
ing mode, we find that C-TDDFT remains quite close to L- 
TDDFT even for large-amplitude oscillations, whereas in the 
breathing mode, C-TDDFT deviates from L-TDDFT by about 
a factor of two for large values of A. 

A more quantitative measure of the deviation of C-TDDFT 
from the exact high-frequency L-TDDFT result is shown 
in Fig. [21 We plot the difference in percent of the 
time average of the absolute value of the scaled power, 
(^2 — ^2)"^ //^ dt\V{t)\/LuA'^, over one cycle for the sloshing 
mode, and over the first and second half-cycle for the breath- 
ing mode, as a function of the amplitude A. As expected, the 
difference is seen to be increasing as the amplitude grows, but 
the two modes exhibit a very different behavior. 

Due to its close kinship to Kohn's mode, it comes as no 
surprise that the deviation is much smaller for the sloshing 
mode, namely at most around 2 % even for very large ampli- 
tudes. The C-TDDFT error of the breathing mode is larger by 
at least an order of magnitude. However, the overall error of 
C-TDDFT in the high-frequency limit compared to L-TDDFT 
remains surprisingly small, as long as the amplitude is not too 
large: For A = 0.2, we get a deviation of about 0.2 % for the 
sloshing mode, and about 5 % for the breathing mode. For 
the largest amplitude considered (A — 0.9), we obtain a de- 
viation of 2.5 % for the sloshing mode, and 20 % and 100 % 
for first and second half-cycles of the breathing mode. In the 
latter case, the deformations are so strong as to invalidate the 



basic assumptions used to derive the simple form of nonlinear 
C-TDDFT. 



2. Non-adiabatic corrections to the ALDA 

The preceding high-frequency analysis shows that the C- 
TDDFT non-adiabatic xc potentials remain close (to within a 
few percent) to the exact L-TDDFT results, except for modes 
with extremely strong deformations such as sloshing modes 
with A > 0.9 or breathing modes with A ^ 0.2. As long 
as the deformations remain within these approximate limits, 
it is reasonable to expect C-TDDFT to be accurate for finite 
frequencies as well. 

To get an impression of the magnitude of the non-adiabatic 
corrections to the ALDA, we plot in figures^jand^Jthe full 
adiabatic and non-adiabatic xc potentials, V^^°^{x,t) and 
V^^°^{x,t) + V^{x,t), where is calculated within C- 
TDDFT using the QV parametrization. Fig. [Olshows results 
for a large-amplitude sloshing mode with A = 0.75, and Fig. 
1 141 for a breathing mode with A — 0.5. The figures show four 
snapshots taken during one cycle of the mode, similar to what 
was presented in Figs. except that now we are plotting 
the total potential. The ALDA xc potential is of course inde- 
pendent of the mode frequency, and we compare it with the 
ALDAh-M potential at three different frequencies: uj ~ O.lWp 
(low-frequency regime), lo — ujp (around cross-over), and 
bj = IQujp (approaching the high-frequency regime). 

The results in Figs. ^]confirm again that the non-adiabatic 
effects for the sloshing mode are relatively modest, even for 
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FIG. 13: (Color online) Snapshots of the adiabatic and non-adiabatic 
xc potentials during one cycle of the sloshing mode, with ampli- 
tude A = 0.75. Black sohd hne: V^^°^. Red long-dashed, green 
medium-dashed, and blue dotted lines: V^^°^ + for frequencies 
u) — O.lwp, Lo — [Dp, and uj = lOwp. 



large-amplitude deformations. The post-ALDA corrections 
become more important for the high-frequency oscillations 
(uj — 10(Dp), where one finds deviations from ALDA of the 
order of 10% at the turning points of the oscillation at r/4 
and 3T/4. For the lower- frequency modes, the non-adiabatic 
corrections to the ALDA stay mostly within about 1%. 

On the other hand, the breathing mode exhibits much more 
dramatic non-adiabatic effects, see Fig. [H] Again we find 
that the post-ALDA corrections are moderate for the lower 
frequencies considered. However, for the high-frequency case 
(io = 10(Dp), we find that at the instances of maximum defor- 
mation (r/4 and 3T/A) the memory effects cause a correction 
to the ALDA of up to a factor of 2, which is indeed substantial. 

The impact of the high-frequency post-ALDA corrections 
is similar for both modes: they tend to oppose the deformation 
of the ALDA potential at the instances of maximum deforma- 
tion of the density distribution. For the breathing mode, this 
means that V^c becomes less broad at T/4 and less deep at 






X (a.u.) 

FIG. 14: (Color online) Same as Fig. 1131 but for the breathing mode 
with amplitude A = 0.5. 



3T/4, and for the sloshing mode, that the potential minimum 
lies closer to the center, and is less sharp. In a more realistic 
calculation where the density, instead of being a given func- 
tion, follows from solving a TDKS calculation with a time- 
dependent driving potential, this would imply that the elas- 
ticity of the electron liquid tends to counteract deformations 
of the density, making the system a more rigid and somewhat 
harder to deform. 

Lastly, to illustrate a case of extreme non-adiabaticity let 
us consider the high-frequency limit of strongly nonlinear dy- 
namics. Fig. 1151 shows a comparison of the ALDA and L- 
TDDFT xc potentials for both modes with amplitude A = 
0.75, for the same snapshots taken during one cycle as pre- 
viously in Figs. E]and^5 For the instants 0, T/2 and T, 
ALDA and L-TDDFT coincide, but for r/4 and 3T/4, large 
deviations occur. The most dramatic non-adiabatic effect is 
observed for the breathing mode at 3r/4, where the elastic 
contribution of L-TDDFT is so large that the resulting xc po- 
tential is of the same magnitude as the ALDA, but with op- 
posite sign. This clearly shows that in situations in which the 
electron density is rapidly and strongly deformed the ALDA 
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FIG. 15: Comparison of ALDA (full lines) and L-TDDFT (dashed 
lines) xc potentials for the sloshing mode (top panel) and the breath- 
ing mode (bottom panel), both with amplitude A — 0.75. The same 
snapshots during one cycle are shown as in figures fT^ and lHI 

becomes drastically wrong, with non-adiabatic corrections of 
the same order of magnitude as the potential itself. 



V. CONCLUSION 

In this paper, we have carried out a detailed comparative 
study of two non-adiabatic local approximations for the xc 
potential in TDDFT, using quasi-one-dimensional model sys- 
tems (whose electronic density is uniform in two spatial di- 
rections, and nonuniform in the third). The goal was to com- 
pare the xc potentials resulting from C-TDDFT and L-TDDFT 
when a given time-dependent density is used as input, and to 
determine the magnitude of the resulting non-adiabatic effects 
when compared to the ALDA. These comparisons were car- 
ried out for two types of collective modes: a sloshing and a 
breathing mode, for a wide dynamical range of amplitudes 
and frequencies. 

In our discussion of non-adiabatic xc effects, we were pay- 
ing particular attention to dynamical regimes where the two 
theories, C-TDDFT and elastic L-TDDFT, are known to be 
exact: C-TDDFT in the linear regime of small-amplitude os- 
cillations, regardless of the frequency, and elastic L-TDDFT 
in the high-frequency limit, for oscillations of an arbitrary am- 
plitude. This helps to shed light on the so far open question 
concerning the validity of C-TDDFT in the nonlinear regime: 
it turns out that, at least in the high-frequency limit, C-TDDFT 



is surprisingly accurate. Our results show that for moderate 
deformations (up to 20% of the initial density distribution) er- 
rors in the xc potential remain within a few percent. This gives 
us good reason to believe that this nonlinear formalism should 
work well for finite frequencies, too (as long as the deforma- 
tions are not too large). 

Let us now address the fundamental question of the mean- 
ing and the importance of "non-adiabatic effects" in TDDFT, 
based on the insights resulting from our study. In analyzing 
these effects, the language of hydrodynamics, as it is used 
in C- and L-TDDFT, proves very useful. In general, non- 
adiabaticity manifests itself through elastic and dissipative 
components of the electron dynamics, which can be distin- 
guished by their characteristic phase lag with respect to the 
ALDA. Both elastic and dissipative effects enter the theory 
via the dependence of xc potential on the deformation tensor. 
Since the deformation is defined relative to the initial state, 
this dependence reflects how much of its history the system 
carries in its memory. A local-in-time dependence on the de- 
formation, which occurs in the high-frequency regime, cor- 
responds to an extremely pronounced memory and a purely 
elastic xc potential. The dissipative contribution formally ap- 
pears in the form of a time-nonlocality in the dependence of 
xc potential on deformations of the electron fluid. In general 
this time-nonlocality tends to shorten the characteristic mem- 
ory time, and, as a result, it somewhat diminishes the elastic 
contribution. In addition it brings about fundamentally new 
effects, such as relaxation and the coiTesponding power ab- 
sorption. Which of the the two contributions (elastic or dissi- 
pative) is dominant, or whether both play an equally important 
role, depends on the dynamical regime in which the system 
under study is evolving. The results from our simple model 
system lead us to the following conclusions: 

Linear regime. For small-amplitude oscillations, dissipa- 
tion is the most important and dramatic consequence of non- 
adiabaticity, which leads to qualitative corrections to the adi- 
abatic dynamics. There is a dynamical range which we call 
"crossover" regime in which the power absorption due to xc 
retardation effects is maximal. This crossover regime occurs 
for frequencies that are comparable to an average plasmon fre- 
quency for the system. On the other hand, for very low or for 
very high frequencies, dissipation vanishes. Elastic effects in 
the linear regime appear to be less important, in a sense that 
they do not qualitatively change an overall behavior of the xc 
potential. However, the quantitative effect of elastic coiTec- 
tions to the linearized ALDA can be significant, especially in 
the high-frequency regime. In particular, for our breathing 
mode the post- ALDA non-adiabatic corrections are in general 
of the same order of magnitude as the dynamic part of the 
ALDA potential itself. 

Nonlinear regime. Here, elastic effects become more im- 
portant, especially at high frequencies, which can lead to sub- 
stantial, and in certain regimes absolutely dominant contribu- 
tions to the time-dependent xc potential. As we discussed for 
the sloshing and breathing modes, the elasticity of the electron 
liquid tends to oppose attempts to subject the electron sys- 
tem to strong and rapid deformations. This general tendency 
also naturally explains why C-TDDFT has been successful for 
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molecular polarizabilities, which are greatly overestimated in 
ALDA.^^ Dynamic polarization of the system corresponds to 
a redistribution of the charge density, i. e., to the deforma- 
tion of an electron subsystem, which causes a counteracting 
xc force. This force is an intrinsically non-adiabatic effect 
that is completely missing in ALDA. 

Our results once again illustrate the special role of Kohn's 
mode in TDDFT. If the electron dynamics sufficiently resem- 
bles Kohn's mode, as is the case for the sloshing mode con- 
sidered here, non-adiabatic effects are generally small. On the 
other hand, for electron dynamics involving high degrees of 
compression, such as our breathing mode, the non-adiabatic 
"corrections" can become several times larger than the ALDA 
potential itself, so that the ALDA completely breaks down, 
leading to a qualitatively wrong behavior A striking illustra- 
tion of such a situation was shown in Fig. 1151 

Finally, most of the discussions in this paper, especially 
those related to dissipative effects, are relevant only for non- 



adiabaticity in extended systems, where purely electronic dis- 
sipation has a well-defined meaning. The situation is less 
clear when one attempts to describe non-adiabatic effects in 
finite systems such as atoms and molecules. Here, the elas- 
ticity of the electron liquid leads to small shifts of excitation 
energies, but one also obtains finite linewidths due to dissipa- 
tion, which is clearly an unwanted effect. The question how 
non-adiabatic xc potentials for small systems should be con- 
structed thus remains an open issue. 
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